Astronomy & Astrophysics manuscript no. 1709 


February 2, 2008 


(DOI: will be inserted by hand later) 





An infrared study of galactic OH/IR stars. I. 
An optical/near-IR atlas of the Arecibo sample*' ** *** 

KM. Jimenez-Esteban 1 ' 2 , L. Agudo-Merida 1 , D. Engels 2 , and P. Garcfa-Lario 3 



in 
o 
o 

(N 

Oi 
<D 

(N 

m 
> 



^1" 

o 

Oh 
6 

03 



13 



1 VILSPA Satellite Tracking Station, Apartado de Correos 50727, E-28080 Madrid, Spain 

2 Hamburger Sternwarte, Gojenbergsweg 112, D-21029 Hamburg, Germany 

3 ISO Data Centre, Science Operations and Data Systems Division, Research and Scientific Support Department of ESA, 
Villafranca del Castillo, Apartado de Correos 50727, E-28080 Madrid, Spain 

Received 22 July 2004 / Accepted 1 1 October 2004 

Abstract. In this paper we present optical and near-infrared finding charts, accurate astrometry (as 1") and single-epoch near- 
infrared photometry for 371 IRAS sources, 96% of those included in the so-called Arecibo sample of OH/IR stars' (Eder et 
al. 119881 Lewis et al, 11990a Chengalur et al.[T99JJ. The main photometric properties of the stars in the sample are presented 
and discussed as well as the problems found during the process of identification of the optical/near-infrared counterparts. In 
addition, we also identify suitable reference stars in each field to be used for differential photometry purposes in the future. 
We find that 39% of the sources (144 in number) have no optical counterpart, 8 of them being invisible even at near infrared 
wavelengths. The relative distribution of sources with and without optical counterpart in the IRAS two-colour diagram and their 
characteristic near infrared colours are interpreted as the consequence of the increasing thickness of their circumstellar shells. 
Among the objects not detected at near infrared wavelengths four non-variable sources are proposed to be heavily obscured 
post-AGB stars which have just very recently left the AGB. Eight additional objects with unusually bright and/or blue near- 
infrared colours are identified as candidate post-AGB stars and/or proto-planetary nebulae. 
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1. Introduction 

Stars with masses between 0.8 and 8 M pass along the 
Asymptotic Giant Branch (AGB) towards the end of their evo- 
lution. Approaching the tip of the AGB they start to pulsate 
and appear as large-amplitude (<; 0.4'" at K-band) long-period 
(> 80 days) variables. The pulsation is accompanied by heavy 
mass loss which forms a circumstellar envelope of gas and 
dust. If the mass loss rate surpasses M > 10~ 6 M year~' the 
dust shell eventually becomes opaque to visible light (Habing 
[19961 . 

Stellar evolution models predict that stars close to the tip of 
the AGB experience thermal pulses (the so-called 'TP-AGB') 
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leading to variations in luminosity and mass loss rates on 
timescales of several ten thousand years. In addition, they may 
change their photospheric chemistry from oxygen- to carbon- 
rich. The determination of the current evolutionary state of a 
particular AGB star is difficult because of the strong variabil- 
ity and the increasing obscuration which limits our observa- 
tional capabilities. Moreover, their pulsation on timescales of 
1-5 years usually requires long observational monitoring pro- 
grams to determine representative parameters. Hence, well de- 
fined AGB star samples, observed over an adequate long period 
of time, are required to test current evolutionary models of this 
short-lived phase of stellar evolution. 

Available samples of AGB stars are in many cases biased. 
Depending on whether the selection is made at optical/near- 
infrared wavelengths or in the far-infrared the samples under 
analysis contain a larger number of low-mass or high-mass 
AGB stars, respectively. Low-mass progenitors are expected to 
experience only moderate mass loss during their AGB lifetimes 
and, thus they will be detected preferentially in the optical and 
near-infrared. An example is the 'dust-enshrouded AGB sam- 
ple' of Jura & Kleinmann (1989 1, which contains stars with 
typical main-sequence masses Mms ~ 1-3 M Q (Olivier et al. 
2001 1. In contrast, AGB stars descending from progenitors with 
larger main-sequence masses (Mms ~ 2 M Q ) probably spend a 
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larger fraction of their TP-AGB life in an obscured state. This 
is a consequence of the higher mass loss rates, which favours 
their detection in the far-infrared. Less biased AGB samples 
covering the complete TP-AGB and all kinds of chemistries 
should therefore include these far-infrared sources, something 
which is possible only since the infrared survey by the IRAS 
satellite was completed (Beichman et al. U988> . 

Samples of oxygen-rich AGB stars ('OH/IR stars') can be 
constructed using surveys of OH maser emission. In such sur- 
veys AGB stars are easily picked up because of their conspic- 
uous double-peaked maser profile. However, these samples are 
also incomplete, as only part of the oxygen-rich AGB stars ex- 
hibit OH maser emission (Lewis & Engels lT995l l. Monitoring 
of early samples (e.g. Baud et al. I198U revealed either very 
long periods or rather weak variations (Engels et al. [1983 
Herman & Habing 11985b . These OH/IR stars were consid- 
ered as massive AGB stars (Mms = 3-8 M on the main se- 
quence), with the long-period variables shortly before and the 
non-variables probably shortly after the departure from the 
AGB. From the analysis of ISO data of galactic AGB and 
post- AGB stars, Garcfa-Lario & Perea-Calderon (2004) have 
recently postulated the existence of three different chemical 
branches of AGB stars, depending on their progenitor mass: (i) 
low-mass oxygen-rich stars; (ii) intermediate-mass carbon-rich 
stars; and (iii) high-mass O-rich stars, where 'hot bottom burn- 
ing' is active at the base of the convective envelope, preventing 
the formation of carbon (instead, the production of nitrogen is 
favoured). It is however not yet settled if low-mass oxygen- 
rich AGB stars evolve into the large mass loss regime observed 
in high-mass oxygen-rich AGB stars, and if so, how long they 
stay in this regime. It is therefore not clear what fraction of 
low-mass stars is present in OH/IR star samples. Actually, the 
relation between AGB star samples selected in different wave- 
length ranges remains to be explored. 

A well-defined sample of far-infrared-selected oxygen-rich 
AGB stars is provided by the Arecibo survey of OH/IR stars 
(Eder et al. fT9l8l Lewis et al. U990al Chengalur et al. fT993t . It 
consists of w 390 IRAS sources, which were detected in the 
1612 MHz OH maser line with the Arecibo radio telescope. 
The sample was obtained from a complete survey of IRAS 
sources with flux densities > 2 Jy at 25 /im, with declination 
0° < 6 < 37° and appropriate colours of AGB stars (Olnon et al. 
1984 1. The OH maser detection qualifies the IRAS source as an 
oxygen-rich AGB star. The positions of the sources in the IRAS 
two-colour diagram are shown in Fig.^ together with a solid 
line which corresponds to the sequence of colours predicted 
for oxygen-rich AGB stars with increasing mass loss (Bedijn 
I1987> that we have named the 'O-rich AGB sequence'. The 
broad range of colours suggests that this sample bridges the 
gap between the optically/near-infrared selected AGB samples 
and the OH/IR stars discovered by OH maser surveys. 

The Arecibo sample is probably composed of objects with 
a variety of variability properties, although the Mira-type large 
amplitude variables are thought to prevail. Unfortunately, vari- 
ability information only exists in the literature for around 50 
sources in the sample, concentrated at the blue end of the IRAS 
two-colour diagram. Most of them are well known Mira vari- 
ables (145 rf < P < 660 rf ) listed as such in the General Catalogue 



of Variable Stars (Kholopov et al. 1998 1 on the basis of opti- 
cal observations. Ten of them have also been monitored in the 
near-infrared (Whitelock et al. 119941 Olivier et al. 2001 ), dis- 
playing amplitudes of 0.5"' < AK < 1.7 m . At least one of these 
blue sources, IRAS 02404+2150, has been reported to show a 
low amplitude semiregular variability (Whitelock et al. 1995 1. 
At the red end a handful of sources has previously been referred 
to in the literature either as variable or non-variable OH/IR 
stars, but usually a detailed knowledge of their light curves is 
still lacking. In addition, a few peculiar sources are also part of 
the sample, such as the very young planetary nebula M 1 -92 
(IRAS 19343+2926). 

In order to study the variability properties of the 'Arecibo 
sample' (periods, amplitudes, and colour variations) we started 
in 1999 a near-infrared monitoring program that we are still 
carrying on. At the end of 2003 we had already completed 1 1 
observing runs with more than 110 nights of observation. A 
subsample is also being monitored optically. The first results 
of this multi-wavelength long-term monitoring program are re- 
ported in this paper. Our ultimate goal is to study in particular 
the oxygen-rich AGB stars with Mms ~ 2 M Q which are proba- 
bly rare among AGB stars discovered optically and in the near- 
infrared, but are common in samples discovered by blind OH 
maser surveys. 

In this paper, which is the first of a series of three, we will 
focus on the determination of the optical/near-infrared coun- 
terparts of all the OH/IR stars in the 'Arecibo sample' and on 
the identification of suitable reference stars in each field to be 
used for differential photometry purposes in the future. The re- 
sult is an atlas of finding charts combining optical images from 
the Digitized Sky Survey (Djorgovski et al. 2001 1 and our own 
near-infrared images, obtained in the J, H, and K bands. 

Sect. 2 describes the observations performed and the acqui- 
sition and data analysis techniques applied. In Sect. 3 we ex- 
plain how we identified the optical/near-infrared counterparts 
and how we determined new accurate coordinates. This is fol- 
lowed by a brief description of the atlas. In Sect. 4 we present 
the near-infrared photometry, while the results obtained are dis- 
cussed in Sect. 5. Finally, the main conclusions are presented 
in Sect. 6. 



2. Observations and data reduction 

Our 'Arecibo sample' of OH/IR stars contains 385 objects, and 
comprises most of the original sources listed in Eder et al. 
il9M , Lewis et al. dl990al l. and Chengalur et al. JT9931 . The 
only objects not considered here were sidelobe responses or 
serendipitous OH maser discoveries, as described in Engels & 
Lewis ( FT996l >. together with three objects (IRAS 18534+0215, 
IRAS 19175+1344, and IRAS 19226+140 1) that were left out 
because they were re-classified by Lewis ( 1997 1 as molecular 
cloud sources. Finally IRAS 19200+2101 was also removed, 
as it turned out to be an erroneous entry in the IRAS PSC (see 
Lewis et al. 20041. 
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Fig. 1. The position in the IRAS two-colour diagram of the stars in the 'Arecibo sample', where the IRAS colours are defined 

F (12} F f25) 

as: [12]-[25] = -2.5 log /ps) an( ^ [25]-[60] = -2.5 log F The continuous line is what we call the 'O-rich AGB sequence' 
(see text). 



2.1. Observations 

The infrared observations under discussion here were per- 
formed during 20 nights in July 1999 and July 2000 at Calar 
Alto Observatory (Almeria, Spain) using the 1.23 m tele- 
scope equipped with the infrared camera MAGIC (Herbst 
et al. 119931 . All sources in the sample, except 13 objects 
not accessible at this time of the year (with A h <a<\2 h ) 
and the very bright, detector saturating, Mira variable RAql 
(IRAS 19039+0809), were observed in the J (1.25 pan), H 
(1.65 /mi) and K (2.20 /mi) near infrared bands. The sources 
observed (N = 371) are listed in Table A2 of the Appendix. All 
except one (IRAS 18549+0208) were observed at least once 
under photometric weather conditions. The MAGIC camera 
uses a 256 x 256 pixels NICMOS3 array that provides an ap- 
proximate field of 5' x5'. Typical on-source integration times 
were 1 minute for the J and H bands, and 30 seconds for the 
K band, resulting in typical upper detection limits of around 
15 m , 14 m and 13 m , respectively. To avoid saturation when ob- 
serving the brightest sources a small aperture configuration was 
used. Standard stars from Elias et al. (1982 1 were observed dur- 
ing each night for photometric calibration. They were followed 
over a range of airmasses to determine the atmospheric extinc- 
tion corrections. 

The main difficulty to observing in the infrared is the high 
sky background level and its fast temporal and spatial varia- 
tion. For the accurate determination of the background we took 
for each object and every filter 5 images with the object cen- 
tered at different positions in the image (dithering or 'moving 
sky' technique). Thus, each measurement consisted of a set of 



5 individual images. Taking the median of the 5 images we ob- 
tained a temporally and spatially well determined sky image 
which was subtracted from the object image. Flat-field images 
were also taken every night for each individual filter in each of 
both aperture configurations to correct for pixel-to-pixel sensi- 
tivity variations during the data reduction process. 

2.2. Data reduction 

During each of the two observing epochs we obtained of the 
order of ten thousand images covering the full sample twice. 
Within the long-term monitoring project we expect to take 
~ 150 000 images, which makes a high level of automation 
in data reduction desirable. Thus, a semi-automated reduc- 
tion procedure was developed combining a self-written IDL 
(Interactive Data Language; Research Systems Inc.) routine 
with pre-existing routines within the IRAF 1 software package. 
The IDL routine first reads in all images taken of in a partic- 
ular observing night from the observing logfile provided by 
the MAGIC camera software. Then it writes an IRAF macro, 
which combines several IRAF routines to perform the data 
reduction which consisted of the following steps. First, a de- 
tector mask was used in order to remove the bad pixels from 
the images. Then these images were divided by the appro- 
priate flat-field images and subsequently the median sky im- 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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ages were subtracted. The sky-subtracted images were then in- 
spected for cosmic ray hits and cleaned by interpolating be- 
tween the neighbouring pixels. Finally, each set of 5 images 
obtained per source and filter were aligned and combined to 
yield the final near-infrared atlas image. 

3. Identification of the optical/NIR counterparts 

In order to determine the optical/near-infrared counterparts of 
the sources in our sample we first considered the best coordi- 
nates available from the literature or from existing catalogues 
and then searched for plausible counterparts at these locations 
both on the optical images from the Digitized Sky Survey and 
on our near-infrared images. 

Improved coordinates with respect to those originally pro- 
vided by the IRAS Point Source Catalogue (accuracy typ- 
ically between 10" -15") were obtained in many cases by 
cross-correlating our sample with the MSX6C Point Source 
Catalogue (Egan et al. 2003), which provides coordinates with 
an accuracy of « 2". Although the MSX survey was limited to 
low Galactic latitudes (< 6° in absolute value), its accuracy is in 
many cases essential for identifying the near-infrared counter- 
part, in particular in crowded regions along the galactic plane 
and for extremely red objects only marginally detectable in the 
K band, as we will see later. 

Out of the 269 objects (73% of the observed sample) lo- 
cated within the region of the sky surveyed by MSX, (essen- 
tially |b| < 6°), we found that 249 objects (93%) had a mid- 
infrared counterpart in the MSX6C Point Source Catalogue 
(see Table \l\ and Table Al in the Appendix). The missing 
objects are all located at the edges of the MSX survey area, 
at galactic latitudes in the range 4.9°<|b|<6°. Two sources 
(IRAS 19161+2343 and IRAS 19206+2517), although not in- 
cluded in the Point Source Catalogue, were identified directly 
on the MSX images, from which we also derived the improved 
coordinates that were later used in the identification process. 

Additional improved coordinates were taken from Lewis et 
al. J1990b> . who determined positions of the OH masers by ra- 
dio interferometry at the VLA for 46 sources in our list. These 
coordinates have errors < 1". Note that part of these objects are 
also MSX detections. 

Once the best coordinates available were determined for 
each source in the sample, we searched for their near-infrared 
counterparts in our images obtained at Calar Alto. For this we 
inspected a 30" x 30" box (or smaller, depending on the ac- 
curacy of the available - IRAS, MSX or VLA - coordinates) 
centered at the nominal Arecibo source position, and searched 
for a plausible (i.e. redder than average) counterpart within the 
box. Usually only one 'red' source was found close to the ex- 
pected position; in many cases this was the brightest source in 
the near-infrared field and/or showed extremely red colours. 

In parallel, we searched also the corresponding position in 
the Second Digitized Sky Survey (DSS2) using the red filter 
image, which covers the spectral range 6 000 - 7 000 A with a 
maximum efficiency around 6 700 A. In many cases it was easy 
to confirm the previous identification made in the near-infrared 
since no optical counterpart (or a very faint one) was seen at 
the corresponding position. 



However, in several cases additional information was 
needed to verify the optical and near-infrared identification. For 
a few Arecibo sources more than one candidate was found with 
appropriate colours. The right counterpart was then uniquely 
identified by searching for variability between two observ- 
ing epochs. The same verification method was used for near- 
infrared counterparts which were detected only in the K band, 
and for which therefore colour information was lacking. 

In most cases (61%), a single, point source counterpart was 
found both in the optical and in the near-infrared. For less than 
one third of the sample (28%), a bright near-infrared counter- 
part was found but so heavily obscured in the optical that noth- 
ing was seen on the DSS2 image above the detection limit of 
S3 20.8'". 11% of the sources are so strongly obscured that in 
the near-infrared the counterpart was not identified above the 
detection limit either (2%), or only in the K-band (9%), being 
completely invisible at shorter wavelengths. 

In all but 8 cases the identification strategy provided a 
plausible counterpart to the Arecibo source, either in the op- 
tical and in the near-infrared or only in the near-infrared. 
The 8 sources without counterpart are listed in Table ^ 
All have MSX and/or VLA coordinates, and their counter- 
parts must have been fainter than our detection limits at the 
two epochs of observation considered. Two of these objects, 
IRAS 18498-0017 (OH 32.8-0.3, a well known long-period 
variable OH/IR s tar, En gels et al. lT983l and IRAS 19440+225 1 
(Lawrence et al. 1990) were reported as detected sources in 
the past but only at wavelengths beyond 3 yum. They do not 
appear to have a counterpart in the 2MASS Point Source 
Catalogue (Cutri et al. 2003) either. Two of the sources in 
Tabled IRAS 18517+0037 and IRAS 19374+1626, were clas- 
sified as variable OH/IR stars in the GLMP catalogue (Garcfa- 
Lario I1992L based on the large value of their IRAS vari- 
ability index (quoted in Column 9). The same classifica- 
tion applies also to IRAS 18475+0353 and the above men- 
tioned IRAS 18498-0017. The remaining three objects have 
a low variability index. One of them, IRAS 18596+0315 
(OH 37.1-0.8, GLMP 862) is a non-variable OH/IR star which 
has already left the AGB (Engels 2002). If the non- or small 
amplitude variability of these objects is confirmed, they are 
probably also non-variable OH/IR stars that have left the AGB 
recently. 

Taking into account that IRAS 18 498-00 17 (OH 32.8-0.3) 
has a period P > 1 500^ (Engels et al. 119831 ) and that the other 
stars in Tableware among the most obscured objects in the 
sample, we conclude that the variable sources are probably 
OH/IR stars with very long periods (> 1 000 rf ) at the very end 
of the AGB, while the non-variable ones must be heavily ob- 
scured post- AGB stars, already in the transition to the planetary 
nebula stage. 

3.1. The atlas of optical and near-infrared 
counterparts 

The selected counterparts are displayed as an atlas of find- 
ing charts. An example of the atlas images is given in Fig. 
|2 For each Arecibo source in our sample a chart was put to- 
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Table 1. Arecibo OH/IR stars without optical/near-infrared counterpart. Coordinates are from radio interferometry (N=2) or from 
MSX. 
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Fig. 2. Atlas images for IRAS 19323+1952. The upper row gives the optical DSS2 R band and the Calar Alto J band images, 
and the lower row the H- and K-band images. The adopted counterpart is marked by a big circle, the reference stars by small 
circles surrounded by squares. 



gether as a mosaic of 4 frames containing the optical image 
taken from the DSS2 in the upper left panel, and the J, H 
and K images from our observations in the upper right, lower 
left and lower right panels, respectively. The size of the field 
shown in each filter is 4.6' x4.6'. For each source we marked 



the position of the optical/near-infrared counterpart with a cir- 
cle in each of the available frames, as well as the position 
of the reference stars used for the astrometric measurements 
(with small circles surrounded by squares). These are the same 
stars which are used for differential photometry purposes in 
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our on-going long-term monitoring program. In those cases 
where an optical/near-infrared counterpart was not found the 
circle was drawn at the position where the source should be 
located according to the best coordinates available (MSX or 
VLA). The complete atlas can be accessed electronically at: 
http://www.edpsciences.org/aa 

3.2. Derivation of improved coordinates 

Once the counterparts were identified, we used the set of refer- 
ence field stars (usually three) visible in both optical and near- 
infrared images to derive new accurate coordinates of the target 
stars. Using the relative distances between the Arecibo source 
and the reference stars in the near-infrared we can determine 
the position of the counterpart with an accuracy equivalent to 
less than one pixel on the optical image. If a source was found 
at that precise position it was considered to be the optical coun- 
terpart of the Arecibo source. 

In a few cases (4%) the near-infrared counterpart was so 
bright that the short integration time used to avoid saturation 
was insufficient to detect any reference star in the field. In these 
cases, a very bright object was also visible at the nominal posi- 
tion in the DSS2 image and this was taken as the optical coun- 
terpart. The plausibility of this identification is supported by 
the blue IRAS colours associated with these extremely bright 
near-infrared sources, which imply that they are not expected 
to be heavily obscured in the optical. Because of the very bright 
near-infrared counterpart of these stars the variability monitor- 
ing is done in the optical to avoid saturation problems during 
their luminosity maximum. 

Taking advantage of the astrometric accuracy (1 "/pixel) 
of the DSS2 images we systematically improved the accuracy 
of the astronomical coordinates of all sources in the sample 
for which either an optical or a near-infrared counterpart was 
found. If the Arecibo source had both an optical and a near- 
infrared counterpart, the coordinates were taken directly from 
the DSS2 image by determining the centroid of the point-like 
emission associated to the optical counterpart. If only a near- 
infrared counterpart was found, the position of the Arecibo 
source was determined from its relative position to the refer- 
ence stars in the near-infrared field, which were always chosen 
to have point-like optical counterparts. In both cases, the typi- 
cal astrometric errors are estimated, from the pixel scale of the 
DSS2 images, to be « 1" both in RA and DEC. The complete 
list (363 objects) of new coordinates is displayed in Table Al, 
where we also list the relative distance from the optical/near- 
infrared counterpart found to the original IRAS coordinates. 
When available, the distance of the associated MSX and VLA 
counterparts to the IRAS position is also provided for compar- 
ison. 

If the Arecibo source had neither an optical nor a near- 
infrared counterpart, then the coordinates assigned are the best 
ones available for the source in the mid-infrared (MSX) or at ra- 
dio wavelengths (for the two with VLA measurements). These 
are listed in Table [l] 



4. Near infrared photometry 

Near infrared photometry was performed on each of the sources 
in the sample for which a near infrared counterpart was found 
considering the whole set of images available. The SEXtractor 
software package (Bertin & Arnouts, 1996) was used in com- 
bination with IDL routines to determine for each source the 
residual background noise cr in each of the five frames avail- 
able per filter and to generate a catalogue with the accurate po- 
sition and the aperture photometry (in integrated counts) of all 
the point sources detected in the field with a signal at least a 3 cr 
above the background. The counts corresponding to the source 
identified as the Arecibo counterpart were then extracted from 
these files and averaged, and a standard error cr' was calculated. 
Individual measurements exceeding the average value by more 
than 3 cr' were not considered for the subsequent flux determi- 
nation and error calculation. The counts were then transformed 
into absolute fluxes taking into account the calibration obser- 
vations made on standard stars which were used to derive the 
zero point of our photometric system and the atmospheric ex- 
tinction. The associated errors were calculated using standard 
error propagation analysis, yielding typical values of ss 0.03"' 
for 'well-detected' sources in the three filters considered. 

In Table A2 of the Appendix we present the near infrared 
magnitudes (or lower magnitude limits) derived for the Arecibo 
OH/IR stars observed. No photometric data are available for 
IRAS 18549+0208, because it could not be observed under 
photometric conditions in any of the observing epochs, nor 
for IRAS 19029+0933 and IRAS 20149+3440, because these 
sources were found to be strongly blended with a field star. Ten 
objects were barely detected in the K band, and were there- 
fore too faint to perform reliable photometry. For them, lower 
magnitude limits in the three filters are given. Because of vari- 
ous technical problems we were not able to derive the K band 
magnitude of IRAS 19422+0933, and the H band magnitude of 
IRAS 20440+0412. We noticed during the data reduction pro- 
cess that several bright OH/IR stars observed with the large 
aperture were measured in the non-linear regime of the detec- 
tor resulting in an underestimation of their brightness. For these 
objects upper magnitude limits are given. 



5. Discussion 

We made a number of consistency checks involving positional 
correlations and the analysis of colours, to validate the near- 
infrared counterparts selected and to single out exceptional ob- 
jects. 

5.1. Positions of the near-infrared counterparts 

We have compared our new positions: i) with those originally 
listed in the IRAS Point Source Catalogue for all the sources 
in the sample with an NIR counterpart; ii) with those for a sub- 
sample of 243 objects with an MSX counterpart; and iii) with 
the accurate VLA radio coordinates of 44 OH/IR stars. In Table 
|2]we list the median and the mean separation in right ascension 
(a) and declination (6) between our Calar Alto coordinates and 
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Table 2. Median and mean separation (and standard deviation) 
between Calar Alto (CA) coordinates and IRAS, MSX and 
VLA ones. 
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Aa (cr A „) 


A6 


AS (cr M ) 






median 


mean 


median 


mean 


CA-IRAS 


363 


3.1" 


5.6" (7.1") 


2.0" 


2.3" (2.2") 


CA-MSX 


243 


1.4" 


1.5" (1.4") 


1.0" 


1.3" (1.4") 


CA-VLA 


44 


1.0" 


1.0" (0.8") 


0.6" 


1.1" (1.4") 
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Fig. 3. Angular separation between Calar Alto and MSX co- 
ordinates for the 243 OH/IR stars in common for which an 
optical/near-infrared counterpart was found. 



the IRAS, MSX and VLA ones, together with the associated 
standard deviations. 

Our new coordinates are in excellent agreement with those 
obtained at the VLA, with an average offset in both RA and 
Dec. of * 1". The only exception is IRAS 18554+0231 with an 
offset of 7" (cf. Table Al). Lewis et al. ( I1990bi list for this 
source a larger astrometric uncertainty than for the others, im- 
plying that the difference is probably due to an error in the radio 
position. 

In addition, there are 243 sources with available MSX co- 
ordinates in our 'Arecibo sample'. The angular separation be- 
tween the MSX coordinates and our new coordinates is shown 
in Fig. [21 The median angular separation is « 1.7". The very 
good positional agreement with the MSX data reinforces the 
reliability of our identifications. 

As expected, much larger deviations are found when we 
compare our new coordinates with the original IRAS coordi- 
nates. Misidentifications are in this case more difficult to de- 
termine, especially if neither MSX nor VLA coordinates are 
available to support our choice of a given counterpart. Fig. 0] 
shows the histograms of these deviations. While in declination 
they are within 6" with just a few exceptions, in right ascen- 
sion we find a considerable number of objects («20%) show- 
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Fig. 4. Absolute coordinate differences in right ascension (up- 
per panel) and declination (lower panel) between Calar Alto 
and IRAS coordinates for the 363 OH/IR stars for which an 
optical/near-infrared counterpart was found. Not included are 
six sources with 30" < \Aa\ < 61". 



ing deviations of the order of 10" and even larger. The larger 
errors in RA is a well known feature of the IRAS Point Source 
Catalogue (Herman et al., 119861 Bowers & Knapp. l 19891 Lewis 
et al., 1990b). The problem is related to the differential accu- 
racy of the IRAS survey in the cross-scan and in-scan directions 
(Beichman et al. [T988l . 

As part of our quality checks we verified that the differ- 
ences found between Calar Alto and IRAS coordinates were 
similar to those found between MSX/VLA and IRAS coordi- 
nates (for the objects with MSX/VLA data available). Large 
but consistent deviations can be attributed to the large errors 
associated to the IRAS coordinates. For the remaining Arecibo 
sources (those with no MSX or VLA counterparts) it is more 
difficult to verify that the coordinate deviations are not due 
to a misidentification. In Table Al there are 14 sources with 
extremely large differences in RA (> 20"). 9 of them have 
MSX coordinates showing very similar deviations with respect 
to the original IRAS coordinates. The remaining five sources 
(IRAS 18033+2229, IRAS 19346+0913, IRAS 20127+2430, 
IRAS 20194+1707 and IRAS 21305+21 18) are bright IRAS 
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as well as near-infrared sources. They also have blue IRAS 
colours and are located at galactic latitudes |b| > 5°, where con- 
fusion is unlikely. We conclude therefore that the identification 
of these counterparts is correct and that deviations with respect 
to the IRAS coordinates up to » 60" are possible, although not 
frequent, and are due to the large errors in the IRAS astrometry. 

5.2. Near-infrared magnitude distribution of the 
'Arecibo sample' 

Fig-0shows the J, H, and K magnitude distribution of the 351 
Arecibo sources with a near-infrared counterpart for which a 
photometric measurement was done at least in the K band. 
Most of the sources observed have bright near-infrared coun- 
terparts with the magnitude distribution peaking around 6.0"' 
in K, 6.5"' in H, and 7.5 m in J. The distributions are however 
not strongly peaked and have a long tail extending to the sen- 
sitivity limit of our observations. We estimate that the distribu- 
tions in Fig.|5]become incomplete in all filters at » 12 m . 18% 
of the observed sources were not detected in the J-band, and 
10% were not detected in the H band. As already mentioned 
in Sect. 13. II for 18 sources (5% of the sample) the counterpart 
was too weak to derive a reliable photometric measurement in 
the K band. 

The observed distribution can in principle be attributed: 
i) to the different intrinsic brightnesses of the sources in the 
sample; ii) to the different apparent luminosities expected from 
sources located in a wide range of distances; and iii) to the dif- 
ferent optical thickness of their Circumstellar Envelope (CSE). 
From the photometric data alone it is unfortunately not possible 
to disentangle the relative contribution of these three effects for 
individual sources. However, considered as a whole, we can 
derive some statistical conclusions from Fig. [5] According to 
previous studies reported in the literature (cf. Habing 1996 1, 
and the results obtained by us (Jimenez-Esteban et al. 2004) 
the distribution of intrinsic luminosities expected in our sam- 
ple of OH/IR stars covers the range between * 2 500 L and 
~26000L Q . This would translate into a dispersion in ob- 
served brightness of the order of « 2.5"'. Moreover, in Jimenez- 
Esteban et al. ( 2004 1 we found that the distribution of luminosi- 
ties is actually strongly peaked around 3 5OO-4OOOL . This 
implies that the dispersion in the observed brightness due to 
the scatter of intrinsic luminosities must be even smaller when 
statistically considered. On the other hand, in Jimenez-Esteban 
et al. (2004) we found that most of our sources are estimated 
to be randomly distributed at distances ranging from 1 to 5 kpc. 
This would mean that their apparent luminosities would cover 
a brightness range equivalent to a; 3.5"'. The combination of 
these two effects is not enough to explain the wide range of 
brightness shown in Fig. [5] (more than 8'" in all filters). Thus, 
the different optical thickness of the CSE of individual OH/IR 
stars must contribute significantly to the observed scatter. 

5.3. The J-H vs. H-K colour-colour diagram 

In Fig. |S]we have plotted the near-infrared colour-colour dia- 
gram J-H vs. H-K of all sources with a detection at least in the 



60 F 
50 ~t 




12 3 4 

H-K 



Fig. 7. The total number of sources in the sample (solid line) 
is compared with the number of sources with an optical coun- 
terpart on the DSS2 images (broken line) as a function of the 
H-K colour. The numbers on top of each bin indicate the per- 
centage of optical detections. 

H and K bands. For objects not detected in J, lower limits for 
the J-H colours were calculated and these are indicated with 
an arrow in this diagram. 

The J-H and H-K colours show an almost linear correla- 
tion extending from (J-H = 0.5; H-K = 0.3) up to (J-H = 4.2; 
H-K = 3.2). The scatter in the range 0.5 < J-H < 3.0 is 0.26"\ 
which is a factor 5-10 larger than expected from the errors 
of our photometry. Beyond J-H = 3.0 the scatter increases to 
* 0.40 m . 

The colour-colour distribution is similar to that obtained 
by Lewis et al. (I2004> from 2MASS data for one third of the 
'Arecibo sample'. At its blue end it encompasses also the re- 
gion where Mira variables are found (Whitelock et al. Q994). 
Unfortunately, the limited sensitivity of our observations espe- 
cially in the J band prevents the full exploration of the red end 
of the distribution. If we extrapolate the correlation found, the 
reddest objects in our sample with H-K > 4 m would have as- 
sociated colours J-H » 6, well beyond the limit of J-H = 4.2 
shown in Fig. [6] 

In general, the position of the sources in Fig.[6]can be ex- 
plained as the result of the combined emission of a cool star 
(T«2500K) and a much cooler dust shell (T<800K). The 
main effect of the shell in the near-infrared colours would be 
to increase the circumstellar reddening, which, together with 
differential interstellar extinction effects from source to source, 
could explain the width of the distribution observed. 

But the correlation shown in the near-infrared colour-colour 
diagram J-H vs. H-K essentially represents a sequence of in- 
creasing optical thickness of the CSEs, where Mira variables 
with still optically thin CSEs are placed in the bluer part, and 
OH/IR stars with thicker CSEs are located in the redder part 
of the diagram. We expect therefore that the fraction of objects 
with optical counterparts will decrease with increasing near- 
infrared colours. 

In Fig. we have plotted the fraction of optically 
identified counterparts on the DSS2 images versus the 
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Fig. 5. J band (dotted line), H band (dashed line), and K band (solid line) magnitude distribution of the sources included in the 
'Arecibo sample' . 



near-infrared colour H-K. As expected, the percentage of 
detections decreases with H-K colour. For H-K % 1 al- 
most all Arecibo sources have optical counterparts, while 
at H-K ^;1.5 the fraction decreases rapidly to <50%. 
The two very blue (H-K < 0.7) Arecibo sources without 
optical counterpart are IRAS 18409+0431 (J = 6.8 m ) and 
IRAS 19178+1206 (H=14.2 m ) (cf. Table Al). The non- 
detection of IRAS 18409+0431 might be caused by variabil- 
ity or extremely high extinction, while IRAS 19178+1206 will 
be later identified in this chapter as a non-variable OH/IR star 
already in the early post-AGB stage. Beyond H-K = 2 the frac- 
tion of optical detections decreases further but does not go 
down completely to zero. The very red stars (H-K > 2.5) with 
optical counterparts are probably highly variable OH/IR stars 
optically detected during a luminosity maximum. Another pos- 
sibility is that the proposed optical counterparts are misiden- 
tified and they are actually faint field stars unrelated with the 
IRAS sources. 

Fig. [S] shows the position of the Arecibo sources in the 
IRAS colour-colour diagram as a function of their near-infrared 
colours. Blue near-infrared sources with H-K < 1.0 are located 
preferentially in the bluer part of the diagram ([12]-[25] < 
-0. 1), while those with a red near-infrared colour (H-K > 2.0), 
most of them heavily obscured sources without an optical coun- 
terpart as we have just shown, occupy the red part of the dia- 
gram ([12]-[25] >0.1). Those with intermediate near-infrared 
colours occupy an intermediate range of [12]-[25] colours, 
as expected. There are however large overlaps between these 
groups. 

The smooth correlation found between near- and mid- 
infrared colours corroborates the generally accepted interpreta- 



tion that both the H-K and the [12]-[25] colour sequences are 
good indicators of the optical thickness of the CSE. The colours 
cover the whole optical thickness range, from the optically thin 
shells of Mira variables to the thick shells of OH/IR stars. As 
the transparency of the shell decreases, stars are expected to be- 
come obscured, initially in the optical range and subsequently 
in the near-infrared at increasing wavelengths. The very red ob- 
jects which could not be detected in the H band represent the 
most extreme AGB population of obscured stars. Most of them 
are actually located among the reddest sources also in the IRAS 
two colour-colour diagram, as we can see in Fig. [8] 

The detection of a few blue near-infrared (H-K < 1 .0) coun- 
terparts with very red [12]-[25] (>0.55) colours could be an 
indicator of the end of the high mass loss regime in these 
stars, which should then be classified as post-AGB stars. Stars 
in the post-AGB phase are expected to reappear again as 
bright sources first in the near-infrared and later in the opti- 
cal range, while the remnant of the circumstellar shell dissi- 
pates. Alternatively, they could also be field stars, with the true 
counterparts being heavily obscured beyond our detection lim- 
its. We will analyse this small group of stars later more in detail. 

5.4. The K-[12]vs. [12]-[25] colour-colour diagram 

Another way to analyse the near- and mid-infrared properties 
of the sources in our sample is to study their distribution in 
the K-[12]vs. [12]-[25] colour-colour diagram (Fig. All 
stars with a reliable K band photometry and good quality IRAS 
fluxes at 12 and 25 //m are plotted. If only upper limits in flux 
(lower limits in magnitudes) are available these are indicated 
with arrows. In stars with optically thin shells in the near- 
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Fig. 6. Near-infrared J-H vs. H-K colour-colour diagram of all Arecibo sources with photometry at least in the H and K bands. 
J-H lower limits are indicated by arrows. The solid line shows the location of black-bodies of different temperatures and the 
dashed lines show the location of the combination of a 2 500 K black-body with dust shells of temperatures 200, 400 and 800 K 
following Whitelock ( 1985 1. In addition we also indicate the reddening vector corresponding to Ay = 10 m . 



infrared the K-[12] colour provides information on the rela- 
tive contribution to the overall spectral energy distribution of 
the near-infrared emission, dominated by the central star and 
by the hot dust surrounding it, and of the mid-infrared emis- 
sion, mainly coming from the cool dust in the circumstellar 
shell. The correlation of the K-[12] with the IRAS [12]-[25] 
colour can then be interpreted as an additional indicator of the 
optical thickness of the CSE of a given source. To calculate the 
K-[12] colour we adopted a zero-magnitude flux in the K band 
of 665 Jy (Koornneef fT983l . 

Although a clear correlation exists between the K-[12] and 
the [12]-[25] colours, as is shown in Fig. [5] there is a very large 
scatter on both axes. Only a small part of this scattering can be 
attributed to the intrinsic errors associated to the photometric 
measurements, mainly in the IRAS data, as the errors associ- 
ated to our K band photometry are in most cases negligible. 
Typical errors are in the 10-15% range at 12 and 25 /mi for 
the sources in our sample which translates into a dispersion of 
just 0.1-0.2'" in the [12]-[25] colour. The influence of inter- 
stellar extinction is also expected to be low, as is indicated by 
the small size of the reddening vector shown in Fig.|9j which 



corresponds to the reddening effect expected for a star located 
at a distance of lOkpc in the K-[12] colour (the effect on the 
[12]-[25] colour would be negligible). The largest contribution 
to the colour dispersion observed is expected to come from the 
strong variability of the sources in our sample. In the case of the 
K-[12] colour, the fact that the near- and mid-infrared obser- 
vations under analysis were non-contemporaneous adds a sub- 
stantial scatter. From data in the literature we know that opti- 
cally bright Miras show typical K band amplitudes of 0.5 - 1 .0'" 
while the most extreme OH/IR stars can reach amplitudes in the 
K ba nd of u p to 4 m (Engels et al. fT983l Le Bertre fT993l Olivier 
et al. 2001 ). In addition, the bluest Miras with OH maser emis- 
sion show a variability at 12 /mi of =s 0.6"' (Smith 2003 1 while 
this variability can reach up to 1 .8"' in the most extreme OH/IR 
stars. This latter value has been estimated by comparing the 
measured flux in 12/im both by IRAS and MSX satellites for 
those Arecibo sources for which these photometric data exist. 
If we include all these uncertainties in the K-[12] vs. [12]-[25] 
colour-colour diagram we obtain a wide area delimited by the 
dashed lines plotted in Fig. [9] where the majority of the stars 
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Fig. 8. IRAS two-colour diagram with the Arecibo sources distinguished according to the H-K colours of their near-infrared 
counterparts. Diamonds (❖) are sources with H-K < 1 .0, crosses (+) are those with 1 .0 < H-K < 2.0, and triangles (A) those with 
2.0 < H-K. Circles (o) correspond to sources with H brightness below the detection limit. The solid line is the 'O-rich AGB 
sequence'. 



are located. This would imply that in most cases the intrinsic 
variability of the sources could explain the scatter observed. 

The K-[12] colour distribution has also been studied by 
Whitelock et al. ( 1994 1 for Mira variables and by Sevenster 
(2002) for OH-selected AGB and post- AGB stars. The Mira 
variables in Whitelock et al. i 1994b cover the blue end of the 
J-H vs. H-K two-colour diagram (Fig. [6} with the majority of 
them concentrated around values of J-H ^1.0 and H-K « 0.5. 
In contrast, the OH-selected AGB and post- AGB stars studied 
by Sevenster ( 2002 1 show a wider distribution in this diagram, 
very similar to what we observe in our 'Arecibo sample' . 

In the K-[12] vs. [12]-[25] two-colour diagram the blue 
oxygen-rich Miras studied by Whitelock et al. ( 1994) are con- 
fined to the region indicated in Fig. |9] In the Arecibo sample, 
K-[12] values are in many cases considerably redder, in agree- 
ment with the values found in the OH-selected sample of AGB 
and post- AGB stars used by Sevenster (2002i. The scarcity of 
sources with very red K-[12] colours in Whitelock et al. (I1994> 
is therefore largely a selection effect. 



For Mira variables the K-[12] colour is interpreted by 
Whitelock et al. 1 1994) as a measure of the mass-loss rate argu- 
ing that the K and 12 /urn fluxes originate respectively from the 
star and the shell. This argument is, however, alleviated with 
the increasing optical thickness of the shell and the correspond- 
ing diminishing contribution from the star. The observed range 
of K-[12] colours for a given [12]-[25] value is therefore a 
result of several contributions and a more thorough study is 
needed to properly explain the observations. This will be pos- 
sible as soon as the K-band variability can be removed from the 
diagram as the result of our on-going monitoring program. 

5.5. Early post-AGB candidates 

Fig. [9] also shows at least four clear outliers display- 
ing a near-infrared counterpart much brighter than ex- 
pected from their red [12]-[25] colours. They are labelled 
with their IRAS name in the K-[12] vs. [12]-[25] dia- 
gram. Namely, viz. IRAS 18520+0533, IRAS 18551+0159, 
IRAS 19343+2926 and IRAS 20160+2734. The most natural 
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Fig. 9. K-[12] vs. [12]-[25] colour-colour diagram, where the K-[12] colour is defined as K-[12] = -2.5 log f^rgf , and the 
[12]-[25] colour is defined as in Fig.^ The dashed lines indicate the dispersion expected from the intrinsic variability of these 
sources. Candidate post-AGB stars are marked with open symbols. The dotted box delimits the region occupied by the type of 
Mira variables studied by Whitelock et al. ( 1994). For details, see text. 



interpretation for their peculiar position in this diagram is that 
they are already in the post-AGB stage. As we have already 
mentioned, once the large mass loss at the end of AGB evo- 
lution ceases, the optical depth of the circumstellar envelope 
decreases and the central star is expected to reappear again at 
shorter wavelengths while the mid-infrared colours redden fur- 
ther for some time (Bedijn 1987, Vassiliadis & Wood 1993 1. 
The cessation of the high mass loss is connected to the end 
of the large amplitude pulsation of regular OH/IR stars. Thus, 
non-variability is expected to be an additional post-AGB stage 
indicator. 

In addition, we have also indicated with open symbols in 
F^ig- El the position of those sources with very red [12]-[25] 
colours and unusually blue near-infrared counterparts that were 
also tentatively classified as candidate post-AGB stars or proto- 
planetary nebulae in Sect. 15.31 As we can see, two of these 
sources, IRAS 18520+0533 and IRAS 20160+2734, are among 
the outliers that we have just identified as candidate post-AGB 
stars, which supports our classification. 

Both samples of candidate post-AGB stars are listed to- 
gether in Table [3] where we give their photometric colours 
H-K, K-[12] and [12]-[25], the IRAS variability index and 
the proposed identification. 

As we can see, several of these stars have already been 
identified as post-AGB stars, proto-planetary nebulae or even 
as planetary nebulae in the literature. 

This is the case for IRAS 18095+2704 a well-known 
high-galactic-latitude post-AGB star, with spectral type F3Ib 
(Hrivnak et al. U988> . for which a detailed spectroscopic anal- 



ysis including determination of chemical abundances can be 
found in Klochkova dl995> . 

The NIR photometry of IRAS 18520+0533 shows no vari- 
ability between our two observing epochs and it has a low 
IRAS variability index. Eder et al. ( 1988 1 found that this source 
has a peculiar 1612 MHz spectrum. Its two strongest peaks 
show the standard profile observed in OH/IR stars, but exte- 
rior to these there is a pair of small flat peaks which may be 
due to a second shell of circumstellar material expanding at a 
larger velocity. This points to the peculiar nature of this OH/IR 
star, but we note that its location in the K-[12] vs. [12]-[25] 
diagram depends strongly on the IRAS 12/^m flux, for which 
there is only an upper limit in the IRAS PSC. 

IRAS 18551+0159 was first identified as an infrared plan- 
etary nebula by Kistiakowsky & Helfand ( 1995) on the basis 
of its strong radio continuum emission at 20 cm and the large 
F([SIII]A9532A)/F(Paschen9^9225A) ratio, which means 
that classification as an ultracompact H II region must be dis- 
carded. The source is extremely red and has no optical coun- 
terpart on the DSS2 image, which is probably a signature of its 
recent departure from the AGB. 

IRAS 19065+0832 (OH 42.6+0.0, GLMP876) is the only 
source in Table [3] showing a large IRAS variability index. It 
was classified in the GLMP catalogue as a variable OH/IR star 
(Garcfa-Lario 1992) on the basis of its large IRAS variability 
index, together with the detection of a strong silicate absorption 
feature in its IRAS Low Resolution Spectrum (it is a 'class 39' 
source). However, this is not consistent with the blue colours 
observed in the near-infrared nor with the non-variability which 
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Table 3. Very red sources with unusually bright and/or blue near-infrared counterparts identified as candidate post-AGB stars or 
proto-planetary nebulae. 
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18095+2704 


0.26 


3.55 


1.11 


15 


Post-AGB star 


18520+0533 


0.93 


3.90 


1.76 


3 


Peculiar OH/IR star 


18551+0159 


1.66 


5.20 


1.87 


25 


Infrared Planetary nebula 


19065+0832 


0.70 


3.17 


1.02 


98 


Non-variable OH/IR star? 


19178+1206 


0.65 


7.89 


0.92 


1 


Non- variable OH/IR star 


19200+1035 


0.82 


1.46 


0.61 





Planetary nebula 


19343+2926 


1.79 


2.15 


1.33 





Bipolar proto-planetary nebula 


20160+2734 


0.29 


-1.04 


0.76 





Semi-regular variable 



is suggested by the very similar photometry in our two epochs 
(only 0.07™ difference in the K-band). Our suggestion is that, if 
the near-infrared source was correctly identified, the star must 
be already in the post-AGB stage, but we cannot discard the 
possibility of a misidentification. This is why our classification 
in Table[3]is only tentative. 

We identify IRAS 19178+1206 (GLMP899) as a non- 
variable OH/IR star in a very early post-AGB stage. The source 
is invisible in the optical while it has a relatively blue near- 
infrared counterpart. The non-variability and the blue near- 
infrared colours are confirmed with our two different epochs. 
However, the near-infrared counterpart of this source might in 
fact be a field star, and the true IRAS source invisible even in 
the K band. 

IRAS 19200+1035 (K3-33)is included as a compact plan- 
etary nebula in the Strasbourg-ESO Catalogue of Galactic 
Planetary Nebulae (Acker et al. 1992 1 although with a wrong 
finding chart. It is clearly identified as a bright source in our 
near-infrared images while it is very faint in the optical. Again 
the non-variability and the blue near-infrared colours are con- 
firmed with our two different observing epochs. 

IRAS 19343+2926 (M 1-92), also known as Minkowski's 
Footprint, is a very well known young, bipolar proto-planetary 
nebula (Bujarrabal et al. U998l l. The photometry for this source 
shows a very low variation between our two observing epochs 
(< 0.04'" in K). 

IRAS 20160+2734 (AU Vul) is classified as a semi-regular 
pulsating star in the SIMBAD database. It was first considered 
a post-AGB star by Jimenez-Esteban et al. (2001 1 during a first 
analysis of our present sample. Lewis et al. (|2004i supported 
this classification on the basis of the erratic detection of its 
different masers, being detected in the 1612 MHz OH maser 
at the sensitivity limit of the telescope by Eder et al. ([£9881, 
non-detected later in a water maser survey by Engels & Lewis 
( 1996 1 or as a mainline OH maser by Lewis (19971), and re- 
cently again detected at 1667 MHz but not at 1665 MHz by 
Lewis et al. i2004t . 

6. Conclusions 

We have presented an atlas of optical/near-infrared finding 
charts and near-infrared photometric observations for 371 ob- 



jects taken from the 'Arecibo sample of OH/IR stars'. Except 
for 8 sources, we successfully identified their near-infrared 
counterparts and determined new positions for each source 
in the sample with an accuracy of xl". The correctness of 
the identifications was carefully assessed using improved posi- 
tional information (MSX, VLA), analysing the near- and mid- 
infrared colour measurements, and in some cases searching for 
variability. The few sources for which no near-infrared counter- 
part was found were identified either as heavily obscured vari- 
able OH/IR stars at the very end of the AGB or as non-variable 
OH/IR stars in a very early post-AGB stage. 

The wide dispersion in J, H and K magnitudes is attributed 
not only to the different apparent brightnesses expected for 
sources located in a wide range of distances and to their dif- 
ferent intrinsic luminosities but also significantly to differences 
in the optical thickness of the circumstellar shells. 

The distribution of the Arecibo sources in the near-infrared 
J-Hvs.H-K colour-colour diagram is interpreted as a se- 
quence of increasing optical thickness of their CSEs, where 
Mira-like variables with still optically thin shells are lie in the 
blue part of this diagram while the more extreme OH/IR stars 
with thicker shells are located in the redder part of the dia- 
gram. Their near-infrared colours can be reproduced with a 
combination of the emission coming from a cool central star 
(T * 2500 K) and from a much cooler dust shell (T < 800 K). 
The dispersion of colours observed along the sequence can be 
explained as differential circumstellar and/or interstellar red- 
dening. 

This interpretation was confirmed by a determination of 
the detection rate of optical counterparts as a function of the 
H-K colour. We found that most of the Arecibo sources with 
H-K< 1.0 are detected in the optical range while those with 
H-K > 2.0 are heavily obscured sources with no optical coun- 
terpart in the DSS2 images. These sources are also located in 
different regions of the IRAS [12]-[25] vs. [25]-[60] colour- 
colour diagram, suggesting that both near- and mid-infrared 
colours are good indicators of the optical thickness of the shell. 

The connection between the near- and the mid-infrared data 
was analyzed with the help of the K-[12] vs. [12]-[25] colour- 
colour diagram. We found a clear correlation between these 
two colours which can also be interpreted as an indication of 
the increasing contribution to the overall spectral energy dis- 
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tribution of the mid-infrared component (emission from the 
cool dust in the circumstellar shell) with respect to the near- 
infrared component (emission from the central star plus the 
hot dust surrounding it) as these stars evolve along the 'O- 
rich AGB sequence'. However, the scatter of this correlation 
is found to be very large, partly because of the strong variabil- 
ity of these sources in the near- and mid-infrared and the non- 
contemporaneity of the observations under comparison. The 
expected variability goes from 0.5 to more than 4 m in the K 
band and from 0.5 to 1.8 m at 12yum, both increasing toward 
redder [12]-[25] colours. Once these uncertainties are consid- 
ered, the distribution observed in the K-[12] vs. [12]-[25] dia- 
gram can be explained as a consequence of the different optical 
thicknesses of the CSEs, with the exception of a few outliers 
which were found to be much brighter than expected in the 
near-infrared according to their [12]-[25] colour. 

The few sources showing very red [12]-[25] colours in 
combination with unusually blue near-infrared counterparts as 
well as the outliers in the K-[12] vs. [12] -[25] diagram have 
been identified as candidate post- AGB stars (some of them are 
well known proto-planetary nebulae or planetary nebulae). 

The proper identification of the IRAS counterparts forms 
the basis for a successful completion of our long-term infrared 
and optical monitoring program begun in 1999, which is aimed 
at determining the variability properties of all the sources in 
the 'Arecibo sample of OH/IR stars' and combine the informa- 
tion collected in the near-infrared with data available at other 
wavelengths in order to study more in detail the photometric 
behaviour of these stars in the context of stellar evolution. 
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